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Abstract 

Unification at Mgut ~ 3 x 10^^ GeV of tlie tliree Standard Model (SM) gauge couplings can 
be achieved by postulating the existence of a pair of vectorlike fermions carrying SM charges 
and masses of order 300 GeV - 1 TeV. The presence of these fermions significantly modifies 
the vacuum stability and perturbativity bounds on the mass of the SM Higgs boson. The 
new vacuum stability bound in this extended SM is estimated to be 117 GeV, to be compared 
with the SM prediction of about 128 GeV. An upper bound of 190 GeV is obtained based on 
perturbativity arguments. The impact on these predictions of type I seesaw physics is also 
discussed. The discovery of a relatively 'light' Higgs boson with mass ~ 117 GeV could signal 
the presence of new vectorlike fermions within reach of the LHC. 
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1 Introduction 



The discovery of the Standard Model (SM) Higgs boson is arguably the single most important 
mission for the LHC Under a somewhat radical assumption that the next energy frontier lies 
at the reduced Planck scale (Mp ~ 2.4 x 10^*^ GeV), it has been found that the SM Higgs boson 
mass lies in the range 128 GeV < rriH < 175 GeV pQ. Here the lower bound of 128 GeV on 
rriH derives from arguments based on the stability of the SM vacuum. More precisely, that the 
Higgs quartic coupling does not become negative at any scale between Mz and Mp. The upper 
bound of 175 GeV or so on mn stems from the requirement that the Higgs quartic coupling 
remains perturbative and does not exceed 47r, say, during its evolution between Mz and Mp. 
Thus, it would appear that discovery of a relatively 'light' Higgs boson (with mass well below 
128 GeV) may signal the presence of physics beyond the SM. 

Supersymmetry is by far the most compelling extension of the SM and its minimal realization 
(MSSM) predicts a relatively 'light' SM-like Higgs boson with mass < 130 GeV. However, in 
the light of LHC, plausible alternatives to supersymmetry deserve careful investigation. For 
instance, it was shown in |2j that the new physics between Mz and Mp associated with type II 
seesaw [3] around TeV scale or higher can yield a 'light' Higgs boson with mass > 114.4 GeV, 
the LEP II bound. A 'light' Higgs boson is also realized in scenarios of gauge-Higgs unification 
with a compactification scale below Mp [1]. 

In this paper we revisit another extension of the SM, proposed several years ago, in which 
new TeV scale vectorlike fermions are introduced in order to implement unification at some 
scale Mqut of the three SM gauge couplings The new vectorlike fermions carry SM gauge 
quantum numbers and their presence therefore modifies the SM Higgs mass bounds based 
on vacuum stability and perturbativity arguments. In particular, by including only a pair of 
vectorlike fermions for which case Mqut — 3 x 10^^ GeV, the vacuum stability bound can 
be lowered from its conventional value of around 128 GeV to a significantly lower value of 
about 117 GeV. To keep the discussion as realistic as possible, we also study the possible 
impact neutrino oscillation physics could have on the Higgs mass predictions. We employ type 
I seesaw for these considerations Note that a more complicated scenario containing several 
new particles (including scalars) can yield Mqut — Mp, with a vacuum stability bound as low 
as 114 GeV. 

This letter is organized as follows. In Section 2, we investigate the impact of new vectorlike 
fermions on the vacuum stability and perturbativity Higgs mass bounds. In Section 3, we discuss 
how neutrino physics utilizing type I seesaw can affect the vacuum stability and perturbativity 
bounds. Our conclusion is outlined in section 4. 
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2 New Fermions and the Higgs Boson Mass 



We introduce the following vectorlike fermions to achieve gauge coupling unification: 

Q (^3, 2, + Q (3, 2, -0 + D (^3, 1, ^) + ^ (3, 1, -0 , (1) 

where the brackets contain the SU (3)c x SU {2)l x U (l)y quantum numbers of the new particles. 
The SM Lagrangian is supplemented by additional terms, and the relevant ones are given by 

Cnew = -KiQD<^' - K2QD<^ - y\Qd^^ - y^q.D^ - yiQu^^' - Mf{QQ + DD) + h.c. (2) 

where $ denotes the SM higgs doublet, = ia2^* its charge conjugate, and we employ the 
standard notation g^, m^, dl for the SM quarks, with i = 1, 2, 3. The parameters 1/123 and ki^2 
are dimensionless couplings. We assume, for simplicity, that the new fermions have a common 
vectorlike mass Mp. As pointed out in [7J, most of the 2/123 couplings have to be very small 
due to constraints from the precision electroweak data. To accommodate this, we will assume 
that the couplings ?/| 2 3 sufficiently small so that they do not give a significant contribution 
in the RGE analysis. However, the ' s allow the new fermions to decay into the SM particles, 
without creating any cosmological problems. 

There are constraints on the ^1^2 couplings and the masses of the new matter fields. The 
most important ones arise from the S and T parameters which severely limit the number of 
additional chiral generations. Consistent with these constraints, one should therefore add new 
matter which is predominantly vectorlike. In the limit where the vectorlike mass Mp is much 
heavier than the chiral mass term (arising from Yukawa coupling to the Higgs doublets), the 
contribution to the T parameter from a single chiral fermion is given by [8j 
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where the Yukawa couplings in Eq. (EJ, f = 246.2 GeV is the vacuum 

expectation value (VEV) of the Higgs field, and counts the number of additional SU(2) 
doublet pairs, which in our case is 3. From the precision electroweak data T < 0.06(0.14) at 
95% CL for = 117 GeV (300 GeV) [9J. We will take < 0.1 as a conservative bound for 
our analysis. We see from Eq. ([3]) that with Mp ~ 500 GeV, the Yukawa couplings Hi can be 
0(1). 

For the SM gauge coupling we employ the two renormalization group equation (RGE) [TU] : 
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where gi {i = 1,2,3) are the SM gauge couphngs and yt is the top Yukawa couphng. 
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For a renormahzation scale fi > Mp, the beta function for gauge couphngs receives an additional 
contribution from the vectorlike fermions, 
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where ^ and Cf^ stand for the contribution which is proportional to the coupling in the 
two loop RGE for gauge couplings. 

For the top Yukawa coupling, we have [10] 
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while the two-loop contribution is given by 
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In solving Eq. (j7]), the initial top Yukawa coupling at /i = Mt is determined from the relation 
between the pole mass and the running Yukawa coupling [TTl [T2], 



6 TV V VT 
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(10) 



with yt{Mt) = \/2mt{Mt)/v and 03 = gl/^n- Here, the second and third terms in parentheses 
correspond to one- and two-loop QCD corrections, respectively, while the fourth term comes 
from the electroweak corrections at one-loop level. The numerical values of the third and fourth 
terms are comparable (their signs are opposite). The electroweak corrections at two-loop level 
and the three-loop QCD corrections are both comparable and of sufficiently small magnitude 
[12] to be safely ignored. 
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For a renormalization scale fi > Mp, according to the Eq. ([2]), the beta function for the top 
Yukawa couphng receives an additional contribution at one loop level as follows: 

6(3^'^ = 3iKl + Kl), (11) 
and the additional two loop contributions are 

) = ^^gl + fgl + 20^3^) {nj + nj) - ^(4 + 4) ' '^^v'ti^ + 4)- (12) 
The one and two loop RGEs for the Yukawa couplings ki and H2 are given by 

-1 { T7^/5i? + TTT^/^i?) • (13) 



Here the one loop contribution is 

= --/i - - 8g' + + 3kI + 3yl (14) 
while the two-loop contribution is given by 

= -^^^9t-^-^9t-l08gl-^^^gl^^^ 

+ (^^7? + fgl + 20^73^) - i2«:i - ^(yf + 4 + yt4 + 44)- 



(15) 



The RGE for the Yukawa coupling K2 is obtained by making the replacement ki O k,2 in Eqs- 
f|T3|) - f|T5|) . This follows from the various quantum numbers listed in Eq. ([1]). As previously 
mentioned, we are neglecting mixing terms involving the new vectorlike particles and the SM 
ones. 

The RGE for the Higgs boson quartic coupling is given by [10] 



with 
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We calculate the Higgs boson pole mass ttih from the running Higgs quartic coupling using the 
one-loop matching condition [T3] . 

According to Eq. ([2]) there are additional contributions to the one and two loop beta 
function for A which are proportional to the ki and H2 couplings. At one loop we have 

6/3'^^^ = 12{kI + kI)X- 12(4 + 4), (19) 

and for two loop 

5/3f = (^Igf - 6Agl^ {4 + 4) - l9li4 + 4) + lOA Q^?? + ^gj + Sgj^ {4 + 4) 

+\9l Q^?? + '^gl^ {4 + 4) - 72A^(«:? + 4) - + 4) + 60(«:? + 4)- (20) 

We next analyze the two loop RGEs numerically and show how the vacuum stability and 
perturbativity bounds on the SM Higgs boson mass are altered in the presence of the new TeV 
scale vectorlike particles. 

We chose the cutoff scale to be Mgut, the scale at which the SM gauge couplings are all 
equal. This choice is motivated by the following argument. Namely, we want to have as much 
as possible model independent analysis and in the realistic GUT's we can have very different 
representation for fields. For instance there are many choice of fileds to break GUT symmetry 
[H], or if one address the question of flavor structure of fermions in the framework of GUT, or 
origin of neutrino mass and etc. Also it is well known that in many GUT the cutoff scale has 
to be very close to the Mgut scale doe to existence of big representation under the GUT gauge 
symmetry, for instance in SO(IO), E(6) etc. 

We define the vacuum stability bound as the lowest Higgs boson mass obtained from the 
running of the Higgs quartic coupling which satisfies the condition A(/i) > 0, for any scale 
between Mz < yU < Mgut- On the other hand, the perturbativity bound is defined as the 
highest Higgs boson mass obtained from the running of the Higgs quartic coupling with the 
condition A(yu) < 47r for any scale between Mz fJ^ Mgut- 

In Figure 1, we present the evolution of the gauge couplings for the SM (left panel) and 
for the extended SM (ESM) containing the vectorlike fermions Q + Q + D + D (right panel). 
As noted in [5], in ESM model with new vectorlike fermions weighing a 100 GeV or so, one 
can realize essentially perfect gauge coupling unification at some scale Mgut- Furthermore, if 
we require gauge coupling unification at a level of around 1% or so, then the new vectorlike 
fermion mass should weigh less than a TeV. For definiteness, we set Mp = 500 GeV in our 
calculation. In this case the SM gauge couplings are unified at Mgut — 3 x 10^^ GeV. As seen 
in Figure 1, the new vector like particles help achieve unification by altering the slopes of the 
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Figure 1 : Gauge coupling evolution in the SM (left panel) and in the extended SM (right panel) . The 
vectorlike mass is set equal to 500 GeV and the gauge coupling unification scale is Mqut — 3 x 10^^ 
GeV. 

three gauge couplings. In particular, the slope of is changed and it becomes larger at Mqut 
in comparison to the SM case. The evolution of the top Yukawa coupling is also affected and 
its value is somewhat smaller at Mqut- 

In Figure 2 we show how the evolution of the two-loop top Yukawa coupling in ESM with 
Mp = 500 GeV. The red dashed line stands for the SM case, and the blue solid line corresponds 
to the ESM with = 0. We also present in Figure 2 the evolution of the Higgs quartic coupling. 
The red dashed line corresponds to the vacuum stability bound for Higgs quartic coupling in 
the SM, and the blue solid line corresponds to the quartic couplings in the ESM. We see that 
at Mqut, the top Yukawa coupling in the ESM is smaller in comparison to the SM case. On 
the other hand, it is well known that in the determination of the SM Higgs boson mass vacuum 
stability bound [Ij, a crucial role is played by the interplay between the top Yukawa coupling 
and Higgs quartic coupling, which have comparable and dominant contributions in the RGE 
for Higgs quartic coupling (See Eq. (fTTll ). The negative sign contribution from the top Yukawa 
coupling makes the Higgs quartic coupling smaller during the evolution. This is how the lower 
bound for Higgs boson mass is obtained in the SM. So, having in the model a smaller value 
at Mgut for the top Yukawa coupling means having a milder contribution in the RGE for the 
Higgs quartic coupling, and this explains why in ESM, somewhat smaller values for the Higgs 
quartic coupling f| can satisfy the vacuum stability bound, compared to the SM. In ESM, the 



lower bound for the SM Higgs boson mass using the one-loop matching condition [13] is found 
to be rriH = 117 GeV, close to the LEP bound of 114.4 GeV [IE]. We estimate a theoretical 



A similar observation was made in ref. |15j when considering the type III seesaw mechanism for neutrinos. 
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Figure 2: Evolution of ttie top Yukawa coupling in the SM (red dashed line) and in the extended SM 
(blue solid line). The evolution of the SM Higgs quartic coupling in the two cases are also displayed. 
We have set Mp = 500 GeV and Ki = 0. 

error in this prediction of about 2 GeV, which is in addition to the errors arising from the 
experimental uncertainties in the determination of the top quark mass and [T7j . 

As mentioned earher, the couphng in Eq. ([2]) can be 0(1) if Mp > 500 GeV. In Figure 3 
we present the Higgs boson mass versus for varying Mp scales. For simplicity, we assume that 
K, = Ki = K,2- The upper solid blue and red curves correspond to the Higgs perturbativity bound, 
and the lower dashed curves correspond to the vacuum stability bound when the vectorlike 
particle mass is taken to be 500 GeV (dashed red) and 1 TeV (dashed blue). It is interesting to 
observe that the perturbativity bound decreases as k increases from zero to k ~ 0.6, and then 
increases as the value of k is increased further. We can easily understand this behavior at one 
loop level. It arises from the interplay between the terms 12A(k^ + k,"^) — 12(k^ + nf) in Eq. 
f|T9|) . Up to K ~ 0.6, the term proportional to k^A dominates over the ~ contribution. So, 
for K, < 0.6, in the RGE in Eq. (fT6|) . we have an effective additional contribution with the same 
sign as the A coupling, which leads to the decrease of the perturbativity bound. For k > 0.6 the 
~ contribution dominates compared to the term k^A, and we have an effective additional 
contribution which has the same sign contribution as the top quark in Eq. (1T6|) . This leads to 
an increasing perturbativity bound as the n coupling increases. Note that we have an upper 
bound K = 0.86 for Mp = 500 GeV, and k, = 0.84 for Mp = 1 TeV. This happens because either 
the top Yukawa or k, coupling becomes nonperturbative before the GUT scale. Corresponding 
to the upper bound for n couplings, we have an upper bound on the Higgs mass: rriH = 191 
GeV if Mp = 500 GeV, and tuh = 189 GeV if = 1 TeV. 
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Figure 3: Perturbativity (solid) and vacuum stability (dashed) bounds on the Higgs boson pole mass 
(mjy) versus k{= ki = K2), with vectorlike particle mass Mp = 500 GeV (red lines) and Mp = 1 TeV 
(blue lines). The maximum value for the perturbativity bound is niH — 191 GeV when k = 0.86. The 
lower bound for the Higgs mass is rriH ~ 117 GeV, with k = and Mp = 500 GeV. 

We see in Figure 3 that the vacuum stability bound gradually increases as the k coupling 
increases. This happens because in the evolution of the Higgs quartic coupling, corresponding 
to the vacuum stability bound, the contribution proportional to the term — k'^ dominates over 
the K^A contribution for lower values of n. So in the RGE for the Higgs quartic coupling (see 
Eq. ( fT6l) ) we have an additional contribution with the same sign as the top quark. This leads 
to the explanation why the vacuum stability bound increases when value of k increases at low 
scale, and they eventually merge with the vacuum stability bound. We obtain the following 
results for the Higgs mass corresponding to the vacuum stability bound: rriH = 117 GeV when 
Mf = 500 GeV, and m// = 119 GeV when Mp = 1 TeV, with k = 0. 

3 Type I Seesaw and the Higgs Boson Mass 

We next consider the impact of type I seesaw physics [6j on the Higgs mass bounds found in 
the previous section. The terms relevant for neutrino oscillations through type I seesaw are 
given by 

C, = -y%kiy]<!>^-^M'i{u'^)Juj + h.c., ^,j = 1,2,3. (21) 

Here U is the lepton doublet, the right handed neutrino, ?/]^ is neutrino Yukawa coupling 
and M]^ denotes the right handed neutrino mass matrix. 
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Figure 4: Perturbativity (solid) and vacuum stability (dashed) bounds on the Higgs boson pole mass 
(mjy) versus k{= ki = K2) in the extended SM, including type I seesaw physics. We consider three 
different type I seesaw scales Mr = 10^^ GeV (red), 10^^ GeV (blue) and 10^^ GeV (green). For 
our calculation we consider a hierarchical neutrino mass spectrum, and we set Mp = 500 GeV. The 
maximum and minimal values for the Higgs mass corresponding to the perturbativity and vacuum 
stability bounds are the same as in Figure 3. 



Above the scale M/j we have the following one loop RGE for = y^, 



dlnfi 16 

The various beta functions are modified as follows: 



/3« ^ /3« + tr [YX] , 
I^S ^ + tr [yX.] , 
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(22) 



(23) 



It is certainly interesting to consider realistic cases of the neutrino mass matrix and mixing 
which reproduce the current neutrino oscillation data. We will consider a scenario in which 
the light neutrinos form a hierarchical mass spectrum. It was shown in Ref. [15] that the 
impact on the SM Higgs boson mass from an inverted-hierarchial neutrino mass spectrum is 
not significantly different from the hierarchial case. 

The light neutrino mass matrix is diagonalized by a mixing matrix Umns such that 

M. = i r^F = UmnsDMLns^ (24) 



2M 
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with Du = diag(mi, m2, ms), where we have assumed, for simphcity, that the Yukawa matrix 
Yj, is real. We further assume that the mixing matrix has the so-called tri-bimaximal form |18] 
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(25) 



which is in very good agreement with the current best fit values of the neutrino oscillation data 
]■ 

For the hierarchical case the diagonal neutrino mass matrix is given by 



^ diag(0, Y^Amfs, \j ^vn^^- (26) 
We fix the input values for the solar and atmospheric neutrino oscillation data as [T9] 

Am^2 = 8.2 X 10"^ eV^ 

Am^3 = 2.4 X 10"^ eV^ (27) 

Our finding are presented in Figure 4 where we plot the vacuum stability (dashed) and 
perturbativity (solid) bound versus k(= ki = ^2), with Mp set equal to 500 GeV. We consider 
three distinct mass scales for the heavy right handed neutrinos, namely, M/j = 10^^ GeV 
(red), 10^"^ GeV (blue) and 10^^ GeV (green). The general picture of the Higgs mass versus k 
coupling is qualitatively the same as in Figure 3. Only the initial values for the Higgs mass 
when K = is taken are changed depending on the type I seesaw scale. According to Eq. fl23l) 
the Dirac neutrino Yukawa coupling Yy gives an additional contribution to the Higgs quartic 
coupling RGE with the same sign as the top quark contribution. It is natural to expect that 
the vacuum stability bound will increase if the Yy coupling is increased. For = 10^^ GeV, 
the vacuum stability bound essentially coincides with the corresponding bounds in Figure 3. 
With Mfj = 10^^ GeV the vacuum stability bound is only slightly altered since Yy is still not 
large at that scale in comparison to the top Yukawa coupling. For M/j = 10^^ GeV we see a 
significant change in the vacuum stability bound since now the coupling Yy is larger than the 
top Yukawa coupling, and the two of them together force the Higgs quartic coupling at low 
scale to be larger in order to satisfy the vacuum stability bound. Note that there is hardly any 
impact of type I seesaw on the perturbativity bound. This is due to the fact that above the 
seesaw scale the Higgs quartic coupling is already larger than Yy. 
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4 Conclusion 



Following ref. [5], we have considered a plausible extension of the SM in which new vectorlike 
fermions carrying SM quantum numbers and with masses of order 300 GeV - 1 TeV are in- 
troduced. This relatively modest extension of the SM, denoted by ESM in the text, leads to a 
rather precise unification of the SM gauge couplings at Mqut ~ 3 x 10^^ GeV, and it also gives 
rise to a vacuum stability bound on the SM Higgs mass of 117 GeV. The perturbativity bound 
on the Higgs mass is estimated to lie close to 190 GeV. The new vectorlike fermions should be 
accessible at the LHC. 
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